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UNCLASSIFIED  ABSTRACT 


A  METHOD  FOR  DETERMINING  PROBABILITY  TR-0158(S3960-30)-2 

DISTRIBUTIONS  FOR  MASS  PROPERTIES  OF  June  1S58 

SYSTEMS,  by  Frank  C.  Bond 

An  analysis  is  performed  to  develop  expressions  for  the  statistical  parameters  of 
system  mass  properties.  It  is  assumed  that  the  component  weights,  centers  of 
gravity,  and  inertias  have  known  probability  distributions.  The  general  forms  of  the 
density  functions  generated  are  deduced  from  the  calculus  of  probability  whenever 
applicable.  Means  and  variances  of  the  system  mass  properties  are  computed  from 
linear  functions  of  the  statistical  parameters  of  the  components.  The  input  random 
variables  may  be  distributed  according  to  any  probability  law  or  combination  of 
probability  laws  for  which  the  moments  are  known.  From  conceptual  design  to 
detailed  specification  drawings,  the  probability  laws  governing  mass  properties 
estimation  change.  This  method  may  be  utilized  at  any  point  in  the  life  cycle  of 
system  design  as  it  is  independent  of  the  particular  forms  of  the  distributions  used 
as  input.  (Unclassified  Report) 
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NOMENCLATURE 


«*>■&> 


NOMENCLATURE  (Continued) 


th 

xi’  ^i’  zi  position  of  center  of  gravity  (eg)  of  i  component  in 

reference  frame  XYZ  (see  Figure  1) 

x,  y,  z  position  vectors  of  system  eg  in  reference  frame  XYZ 

$  angle  between  the  X  axis  and  that  principal  axis  having 

the  least  moment  of  inertia 

$  ,  Q  projected  angles  of  principal  axis  offset  angle  from  the 

xy  xz 

reference  axis,  X,  in  the  XY  and  XZ  planes  respectively 
estimate  of  fi  from  simulation 
estimate  of  variance  from  simulation 
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SECTION  ! 

INTRODUCTION 

The  mass  properties  of  a  system  are  functions  of  the  mass  properties  of  the 
constituents  comprising  that  system.  If  the  component  properties  are  known,  then 
the  system  properties  are  uniquely  determined.  Consider  a  system  having  components 
whose  mass  properties  are  random  variables.  If  these  random  variables  are  com¬ 
pletely  defined,  then  the  system  mass  properties  are  random  variables  having  unique 
probability  laws.  The  problem  here  is  to  determine  these  system  mass  properties 
probability  laws  given  the  statistical  distributions  of  the  components. 

In  this  paper,  a  model  describing  the  statistical  nature  of  mass  properties  of 
systems  is  developed.  Given  the  probability  laws  governing  the  weights,  centers  of 
gravity  (eg),  and  inertias  of  the  individual  components,  the  means  and  variances  cf 
the  system  mass  properties  may  be  computed.  Complete  information  about  these  laws 
is  not  required;  only  the  moments  of  the  input  random  variables  are  utilized  as  input. 
Appendix  A  lists  some  of  the  well-known  probability  laws  with  their  associated 
distributions  and  moments. 

A  Monte  Carlo  simulation  (Appendix  B)  was  used  to  verify  the  theoretically 
derived  expressions  for  mean  and  variance  for  all  system  mass  properties  of  a 
typical  reentry  vehicle.  The  component  mass  properties  were  all  input  as  uniform 
random  variables  for  convenience.  The  theoretical  results  were  well  within  the  errors 
associated  with  estimating  parameters  from  finite  samples  (a  sample  size  of  100 
was  used).  Furthermore,  all  mass  properties  appeared  to  be  approximately  normally 
distributed  regardless  of  whether  the  Central  Limit  Theorem  could  be  applied. 
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SECTION  !! 

THE  MASS  PROPERTIES  MODEL 


Consider  a  physical  system  having  n  components.  The  weight  of  the  i 


.th 


component  is  w..  The  center  of  gravity  of  the  i 


.th 


f 


component  is  located  at 
(x.,  y.,  z.  \  within  the  system  reference  frame  X  Y  Z  (see  Figure  1).  Let 


X.  Y.  Z.  be  a  rectangular  coordinate  system  whose  axes  are  parallel  to  X  Y  Z 
i  i  i  ^ 

and  whose  origin  is  at  the  center  of  gravity  of  the  i  component.  The-moments 

and  products  of  inertia  of  w.  with  respect  to  X.  Y.  Z.  are  IQX  ,  I  ,  I  and 

i  oyi  °  i 


I  ,  I  ,  I  ,  respectively, 
oxy.’  oxz.5  oyz.5  r  J 


i  i  •'l 

Using  these  component  mass  properties,  the  system  mass  properties  with  respect 
to  the  arbitrary  reference  frame,  XYZ,  may  be  computed  as  follows: 

Weight 


w 


=  £ 

i=l 


w. 

l 


Center  of  Gravity,  x  Coordinate 


n 


X)  w.  x. 

-  fe  11 


x  = 


n 

E  w 

i=l  1 


,  similarly  for  y  and  z 


Moment  of  Inertia  about  X 


ft  /  2  2  \  ^ 

I  -  Y\  w.  ( y.  +  z.  )  +  V  I  .  similarly  for  I  and  I 
xx  *-i  i\Ji  l  /  ox.  ’  J  yy  zz 

i=l  i=l  i 


3 


2*  0718 


similarly  for  d 


Principal  Axis  Offset  Angle 


6  = 


tan 


6  +  tan2 


X  PRINCIPAL  AXIS 


Figure  2.  Principal  Axis  Offset  Angle 
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SECTION  III 

PROBLEM  STATEMENT 


Determine  the  probability  density  functions  of  the  following  system  mass 
properties: 


Weight 

w 

Components  of  Center  of  Gravity 

x,  y,  z 

Moments  of  Inertia 

*xx’  *yy’ 

Products  of  Inertia 

I  ,  I 
xy’  xz 

Components  of  Principal  Axis  Offset  Angle 

®xy’  ^xz 

Principal  Axis  Offset  Angle 

e 

In  particular,  for  each  of  the  above  random  variables,  determine  the  mean  (//), 
variance  (<x2),  and  the  general  form  of  the  probability  densiiy  function  (normal, 
uniform,  etc.)* 
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SECTION  IV 

•  i 

RESULTS 

A.  GENERAL  FORMS  OF  THE  SYSTEM  MASS  PROPERTIES 
DISTRIBUTION  FUNCTIONS 

For  sufficiently  large  n,  the  Central  Limit  Theorem  states  that  the  weight  and 
the  moments  and  products  of  inertia  are  approximately  normally  distributed.  Since  the 
equations  for  the  center  of  gravity  and  principal  axis  offset  involve  dependent  terms* 
the  theorem  cannot  be  applied  to  these  latter  mass  properties.  However,  the  simulation 
(Appendix  B)  verifies  that  the  dependency  of  these  terms  is  weak  enough  for  them  to 
still  be  normally  distributed  for  most  practical  purposes.  A  good  rule  to  follow  when 
no  other  information  is  available  is  to  assume  that  all  sums  of  random  variables  are 
normally  distributed.  If  a  more  precise  density  function  is  required  in  order  to  make 
probability  statements,  simulation  may  be  used, 

B.  THEORETICALLY  DERIVED  EXPRESSIONS  FOR 
MEANS  (fi)  AND  VARIANCES  (O2) 

The  following  results  are  in  terms  of  the  first  four  moments  ^E  [•] ,  E[* 2], 

E[-3],  E  [•  4])of  the  component  mass  properties.  These  means  and  variances  are 
independent  of  the  forms  of  the  system  distribution  functions. 

Weight,  w 

"w  =!  Eh] 


X  Coordinate  of  Center  of  Gravity,  x 


where 


=  £  (Eh2]Eh2]-E2WE2[\ 


and 


Cov  [t,  w]  =  £  E  S  [wj  E  [x.]  E  [w.] 


£  E[w.2]  e[x.]  -„w  (f  E  [w.]  E  [x.]j 


Replacing  x^  by  y.  and  z.,  we  get  the  expressions  for  the  means  and  variances 
of  y  and  z  respectively. 
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Replacing  y.  by  z.  in  the  above  equations  gives  the  parameters  for  I 

*  *  3 


XZ 


Angle  of  Principal  Axis  Offset  Angle  in  XY  Plane,  e 


2L 


12 


Replacing  by  z.  in  the  above  equations  yields  the  expressions  for 
Principal  Axis  Offset  Angle,  6 


e  . 

xz 


' e 


+  a , 


xy 


G 

xz 


e 


+  (*, 


xy 


e 


xz 


where 
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SECTION  V 
THE  ANALYSIS 


A.  DISTRIBUTION  OF  w 

By  the  Central  Limit  Theorem,  for  independent  w.  and  large  enough  n, 
we  know  that  w  is  normally  distributed  and  has  a  mean 


It  has  been  shown  that  w  converges  to  a  normal  random  variable  rapidly 
even  for  n  as  small  as  3  or  4  when  the  E  jw.j  are  of  the  same  order  of  magni¬ 
tude.  A  mass  properties  probability  distribution  simulation  (Appendix  B)  of  an 
operational  reentry  vehicle  with  n=7  showed  the  weight  to  be  approximately 
normally  distributed  even  though  the  E  |w.J  were  not  all  of  the  same  order  of 
magnitude.  In  the  simulation,  the  w.  were  assumed  uniformly  distributed, 
which  represents  a  conservative  test  for  w  approaching  normality. 
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B.  DISTRIBUTION  OF  x 


Let 


t. 

x 


W.  X. 
1  X 


Since  w.  and  x.  are  independent, 


and 


E  [‘i]  -  E  [wi]  E  [x>] 

eM=  eM  eH 


Using  the  identity 


y  M =  E  M  ■  e2  h 


we  get 


V  |ti 


[‘i] =  E  [wi2] E  [xl2]  -  e2  [wi]  e2  [xi] 


Now  let 


n  n 

1  -  ?  tj  -  Z  W  x 

1=1  1=1 


Since  each  of  the  t.  are  independent,  we  have  the  mean  of  t 

\  -  4]  ■  E[s  =  |  4] 
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or 


and  the  variance  of  t 


Expanding  x  as  a  function  of  t  and  w  in  a  Taylor  series  about  fi^  and 
Mw  respectively, 


-  t  _  "t  *  -  "t 
x-w-W 


"w 


) 


+. . . 


Using  only  the  first  three  terms  of  the  expansion,  the  mean  of  x  is  found  to  be 


Combining  the  last  two  expressions  and  squaring, 


-2 

x 


2  ft 

=  "xJ+  - 

-jL 
..  2 

H-2 

*  ~2  (w  '  "w)' 

u 

'w 

//  / 
'w 

2(1-  . 

-  ^t) 

2/i-- 

+  -  S  [ t 

—  -  It  _  II  \  / 11/ 

+  n  V 

2  (W 
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The  reentry  vehicle  mass  properties  simulation  (Appendix  B)  showed  x  to  be 


approximately  normal  despite  the  dependency  of  the  terms 


xiwi 

w 


.  The  theoretical 


expressions  above  for  and  a—  agreed  very  well  with  the  sample  means  and 

X  X 

variances.  Using  y.,  z.  in  the  above  equations,  we  get  the  equations  of  the  means 
and  variances  of  y  and  7  respectively. 


C.  DISTRIBUTION  OF  I 


xx 


The  mean  of  I  can  be  found  immediately  as 


Ur 


=  E 


XX 


['»]  ■  ‘ft  [’•’  *  •■■)  *  £ 

•S-H  (*M •*['*])•  £*n 


Assuming  that  the  terms  are  all  independent,  we  have  the  variance 


—  V 
xx 


W  ■  fi 


(  2  2 
Wi  \yi  +  zi  } 


•  s  VN 


where 


v[w.  (y.2  +  z.2j]=  E  [wj2  (y.2  +  z,2 )  ]  -EZ[w.  (y/  +  z.2)] 

=  eM(eM  +  2EH  eM+ 
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Hence, 


n 


xx 


-  E 


i=l 

2 


!h2]  (eM +  2E[yi2]  Eh2] +  eM) 

h  (eM + eM)2 


Note  that  and  a  ^ 

xx  XX 


are  exact.  The  Central  Limit  Theorem  would 


/  2  2\ 

guarantee  that  I  is  asymptotically  normal  if  the  terms  w.  I  y.  +  z.  )  and 

I  were  independent.  It  was  shown  in  the  previously  mentioned  simulation  that 
x 

I  is,  for  all  practical  purposes,  approximately  normal.  The  statistical  parameters 

for  I  and  I  are  found  similarly, 
xy  zz  J 


D.  DISTRIBUTION  OF  I 


xy 


The  mean  of  I  is 


%  ■  “!v!  ■  ‘If, 

■  S  (*H  *H  *W  *  "[w,]) 


Assuming  each  term  independent,  we  find  the  variance  of  I  to  be 


- 

?  =  v[l  1  =  V  T  /w.  x.  y.  +  I  \ 

!xy  l  xyj  iti  l  1  1  1  oxv 


XI  r  -J 

£  V  w.  x.  y.  +  I  I 

[ 1  1  1  °^i\ 
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p  i. 


or 


and  the  variance  of  D  is 

xy 


7  2  =  V  [d  1  =  V  T)  w.  ( x.2  -  y.2)  +  I  -  I  I 

Dxy  l  **1  '  [fi  |  1  '  1  '  %  °xij 


Assuming  the  three  terms  are  independent, 


V%S(vh^2-y^]  +  vN-vN) 

or  , 

V =  £  fM  (eM  •  »M  eM + eM) 

-e2w(eM-eM2 

+  e[i  2]  -  E2[l  1  -  E'fl  21  +  E2fl  I 

L  oyi  j  [oyij  Loxi  J  [oxiJ 

Expanding  6  as  a  function  of  I  and  D  in  a  Taylor  series  about 
xy  xy  xy 

the  points  1^  and  respectively, 


I  7 1 

0  =  _SL  =  — SSL.  + 

>  Dxy  *D 

xy 


The  mean  of  0  is 
xy 


'i  -  "i 


-  a  r-  ■  %  +- 


°xy  '  V 
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Assuming  the  component  moments  and  products  of  inertia  are  independent,  we  have 


Replacing  y  by  z  in  the  above  equations  yields  the  parameters  for  9  . 

xz 


F.  DISTRIBUTION  OF  6 

For  small  values  of  S,  we  can  rewrite 


24 


I  I 
xy  xz 


D  D 

xy  xz 


We  can  then  write 


e  e 

xy  xz 


is- Via 

,Dxy)  \Dxz 


The  first-order  approximation  for  the  mean  of  0  6^  is 


E  \d  dl=i 
xy  xz  UF 

l  J  yz 


where 


=  T]  y]  E  I  w.l  Elx.]  E  y.  1  E  w.  j  Ex.  E  z.[ 

m  j=i  H  H  L\l  uJ  L 3j  bJ 

*  s  eM  42]  eH  4J 
+(s  eNK  +(s  eMK 


SeH)(s  eN) 


’  V 


and 


\ ‘ k..t  eH (*M  - eM) Eh](Eh1  - eM) 

i2]  (Ek) + Ekl)+  eH  e[ 


n 

+  E  E 

i=l 


x.  -  E  x 
1 


+  £  I  E 

xy  i=l 


+  "d  E* 

XZ  1=1 


n 

-fE  * 

ii=l 


n 

* 

I 

oz. 

1 

-  E 

i=l 

E 

I 

ox. 

i 

n 

r 

I 

oyiJ 

-  E 

i=l 

E 

I 

OX. 

L  iJ 

\ 


* 

n 

* 

\  /  n 

r  i 

n 

■  " 

s 

-  E  e 

i=l 

*OX. 

l 

SE 

•r-*' 

-Ee 

i=l 

*ox. 

l 

The  covariance  term  can  then  he  computed  as 
CovT#  ,  e  1  =  E  Tfl  0  1  -  u  H 

L  xy’  xzj  |_  xy  xzj  i 


9 

xy  xz 


In  some  cases  the  computation  of  Ha  may  be  considerably  simplified  by 

U 


noting  that 

Cov[v  eA  - '%  ^ 


where  p  is  the  correlation  coefficient  between  0  and  0  ,  and  is  bounded, 

xy  xz’  ’ 


P  < 
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If  it  happens  that 


2  2  „2  2 

17  xz  ^xz  xy  xz 


»2 


"0  "0 

xy 


°0  a0 
xz  xy  xz 


(i.  e. ,  twice  the  magnitude  of  the  second  term  is  negligible  compared  to  the  magnitude 
of  the  first  term),  then  we  can  write 


2 

ae 

2  2 
+  "0  *0 

>> 

2fif,e  "«v  ''a  ff« 

xz 

xy 

xy  xz 

pxy  xz  xy  xz 

or 


^0 

2 

*0 

al 

>> 

2  V  "e  Cov  [v  »«] 

xz 

xy 

xy 

xz 

xy  xz  L  J  J 

and  Hq  is  written  neglecting  the  covariance  term 


U  s  lu%  +  +  — — 

o  V  v  / 


2  2  .  „2  2 

ae  ** e  °o 

xz  xy _ xy  xz 


2  2 

e  M  0 

xy  xz 


3/2 


The  variance  for  9  may  be  readily  found  since  we  have 


M- 


2  +  2 
xy  xy 


and 


E 


e  2  +  u«  2 

XZ  XZ 
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The  variance  is 


or 


°e  2  +  °e  2  +  "e  2  +  Pe  2 

xy  xz  xy  xz 


2 


The  theoretical  parameters  for  6  agreed  quite  well  with  those  found  by  simulation 
for  the  reentry  vehicle.  The  distribution  function  was  shown  to  be  very  nearly  normal. 
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SECTION  VI 


CONCLUSION 

Any  sum  of  independent,  normal,  random  variables  is  normally  distributed. 

For  large  n  ,  the  Central  Limit  Theorem  states  that  the  sum  is  approximately 
normal  regardless  of  the  forms  of  the  independent  component  distributions.  The 
largeness  of  n  must  be  measured  primarily  as  a  function  of  those  components  which 
dominate,  deviation  from  the  theoretical  normal  will  also  occur  if  the  terms  are 
dependent.  If  the  normal  model  is  inadequate  for  a  particular  problem,  a  simulation 
may  be  performed  to  establish  the  general  form  of  the  distribution  function.  In  the 
absence  of  other  information,  the  simulation  shows  that  the  normal  distribution  is  a 
reasonably  good  first  approximation. 

The  derived  means  and  variances  were  verified  (Appendix  B).  Although  some 
of  the  equations  seem  formidable,  they  involve  only  simple  algebraic  manipulations 
and  therefore  are  easily  programmed  for  computer  analysis.  The  equations  yield 
good  approximations  of  the  true  means  and  variances  and  are  valid  for  all  possible 
component  and  system  distributions. 

The  form  of  the  probability  density  function  used  in  conjunction  with  the  derived 
mean  and  variance  completely  describes  the  statistical  nature  of  the  mass  property  in 
question.  Care  must  be  used  in  making  probability  statements,  however,  because  the 
density  function,  the  derived  parameters,  and  especially  the  input  distributions  are 
really  estimates  of  reality.  For  example,  the  simulation  used  components  which  were 
uniformly  distributed.  It  would  be  impossible  to  make  a  3a  probability  statement 
about  the  weight  based  on  the  theoretical  normal  because  +  3aw  exceeds  the 
largest  possible  value  w  may  have.  The  reason  for  this  is  that  the  input  components 
are  uniform  and  therefore  have  finite  upper  bounds.  If  we  had  assumed,  instead,  that 
the  components  were  normal,  then  the  3a  statement  would  have  been  exact. 
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APPENDIX  B 

MASS  PROPERTIES  SIMULATION 


In  order  to  verify  the  theoretically  derived  expressions  for  mean  and  variance 
and  to  gain  some  insight  into  the  general  forms  of  the  probability  density  functions,  a 
simulation  was  performed.  The  object  considered  is  a  typical  reentry  vehicle  sub¬ 
divided  into  seven  components.  Because  the  simulation  was  performed  by  hand,  the 
component  mass  properties  were  assumed  to  be  uniformly  distributed.  For  each 
component  mass  property,  a  sample  of  size  100  was  generated  according  to  its 
prescribed  probability  law.  The  system  equations  were  then  used  to  generate  samples 
of  size  100  for  each  system  mass  property.  Following  this,  estimates  for  the  means 
and  variances  were  derived  and  compared  to  the  theoretical  results  of  the  same  model. 
Minimum  variance  unbiased  estimates  were  used.  For  the  mean  and  variance,  they 
are  of  the  form 


A 

u 


rnj=l 


X. 

J 


and 


m 


^  2 

a 


m 


respectively.  Here  x^  represents  the  j 


•th 


sample  of  the  arbitrary  mass  property,  x. 
The  sample  size  is  m.  In  order  to  display  the  verification  of  the  theoretical 
parameters  without  releasing  classified  information,  the  following  criterion  was  used. 
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A  /\  2 

The  estimates  u  and  a  are  random  variables  with  approximate  standard 
deviations  as  follows: 


^2 

The  approximation  results  from  the  fact  that  the  point  estimate  of  cr  is  used  in 
place  of  the  unknown  variance.  The  expressions 


a2 

SDCcr2] 


and 


SD[u] 


give  a  measure  of  the  difference  between  theory  and  simulation  in  terms  of  the 
standard  deviations  (inherent  errors)  of  the  estimates.  These  results  for  the 
simulation  are  listed  in  Figure  B-l. 

Distribution  functions  resulting  from  the  simulation  are  plotted  in  Figures  B-2 
through  B-7.  They  are  transformed  to  zero  mean  and  unit  variance  for  security 
reasons.  Theoretical  normal  curves  with  zero  mean  and  unit  variance  are  super¬ 
imposed  for  comparison. 
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Figure  B-l .  Difference  befween  Theoretical  and  Simulated  Parameters  in 
Terms  of  Standard  Deviations  of  the  Estimators 
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